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Abstract. Single electron transistors (SETs) are very sensitive electrome- 
ters and they can be used in a range of applications. In this paper we give 
an introduction to the SET and present a full quantum mechanical calcu- 
lation of how noise is generated in the SET over the full frequency range, 
including a new formula for the quantum current noise. The calculation 
agrees well with the shot noise result in the low frequency limit, and with 
the Nyquist noise in the high frequency limit. We discuss how the SET 
and in particular the radio-frequency SET can be used to read out charge 
based qubits such as the single Cooper pair box. We also discuss the back- 
action which the SET will have on the qubit. The back action is determined 
by the spectral power of voltage fluctuations on the SET island. We will 
mainly treat the normal state SET but many of the results are also valid 
for superconducting SETs. 
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1. Introduction 

The single-electron transistor (SET) is known as a highly sensitive elec- 
trometer [1, 2] based on the Coulomb blockade [3, 4]. Electrons tunnel 
one by one through two small-capacitance tunnel junctions, with capaci- 
tances Cl and Cr. The gate, with capacitance C g , is used to modulate 
the generated current, and the island is also capacitively coupled to the 
system to be measured through C c . The charging energy Eq = e 2 /2Cs, 
Cy, = Cl + Cr + C 9 + C c , associated with a single electron prevents sequen- 
tial tunneling through the island at voltages below a threshold Vt, which 
can be controlled by applying a voltage V g to the gate. This Coulomb block- 
ade is effective at temperatures T < Ec/ks and for junction resistances 
larger than the resistance quantum Rk = h/e 2 ~ 25.8kf2. 

SETs can be made using several different technologies, for example 
they can be fabricated from metallic (often aluminum based) systems, from 
quantum dots in GaAs and silicon, or from carbon nanotubes. 

Often the SETs are operated at low temperatures, however room tem- 
perature operation has been demonstrated in several cases [5, 6, 7]. 

Both the ultimate sensitivity and the back action of SET during mea- 
surement is determined by the noise in the SET. By understanding the 
noise in the SET we can optimize the use of the SET for each application. 

1.1. RF-SET 

The conventional SET, based on measuring either the current or voltage 
across the transistor, has suffered from the relatively large output resistance 
R of the transistor. For the typical resistance values of 100 kQ and cable 
capacitance of C ~ InF, the corresponding RC time limits the bandwidth 
to a few kHz. 

With the invention of the radio-frequency SET (RF-SET) [8] the SET 
was made fast and very sensitive. By connecting the SET very close to a cold 
amplifier the upper frequency limit was improved to about 1 MHz[9, 10]. 
With the invention of the RF-SET [8] frequencies above 100 MHz could be 
reached. 

The operation principle of the RF-SET is similar to that of a radar. A 
weak RF-signal (the carrier) is launched via a directional coupler and a bias 
tee, towards the tank circuit in which the SET is embedded. The reflected 
signal depends critically on the dissipation in the tankcircuit, and thus in 
the SET. This means that the carrier is modulated by the gate signal. 
The reflected signal is amplified by a cold amplifier and a number warm 
amplifiers. The reflected signal can be analyzed either in the frequency 
domain or in the time domain. Typical carrier frequencies are in the range 
of 0.3-2 GHz. 
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Figure 1. 
(SCB). 



Schematic figure of the SET capacitively coupled to a Single Cooper-pair Box 



Charge sensitivities of 1.2 • 10 5 and 3.2 • 10 6 e/VWz have been reached 
for signal frequencies of 100 and 2 MHz, respectively [8, 11]. The RF-SET 
can be used as a readout device in applications from very sensitive charge 
meters and current standards [4] in which electrons are counted or pumped 
one by one, to read out of quantum bits[12, 13, 14, 15] or to work as photon 
detectors [16]. 



2. SET - orthodox theory 

Consider a small metallic SET island coupled via low transparency tunnel 
barriers to two external leads, and coupled capacitively to an object to be 
measured. In Fig. 1 we have used the example of a Single Cooper-pair Box 
(SCB) [17] which may serve as a qubit in a quantum computer. The SCB is 
controlled by a voltage source V£ . 

Following orthodox SET theory[2, 3] we use the integer number of elec- 
trons on the SET island (N) as a basis for describing its dynamics. This is 
motivated by the low transparency junctions. The charge on the island is 
the sum of the electrons, the background charge, and the charge induced 
by the the voltages on the three island capacitances. Electrostatics gives 
for the Coulomb charging energy E c h and island potential Vj: 
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where n x e = V^*C| 6 is the charge induced on the gate capacitance C* b 
by the external voltage source V£ b . a L / R describes the capacitive voltage 
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division over the L/R junction, were we have neglected the small gate and 
coupling capacitances {C c ,C g } <C {Cl, Cr}. 

The orthodox theory neglects mainly three things: the effect of the 
electromagnetic environment of the SET [18], higher order tunneling pro- 
cesses, i.e. cotunneling [19], and that for high frequency dynamics the 
transition rates are frequency dependent. As long as the electromagnetic 
environment has a low impedance compared to the quantum resistance 
Rq = h/4e 2 6.45kS7 (which is often the case for SETs) the corrections 
due to the environment are small. The same is true for the cotunneling cor- 
rections as long as the resistance of the SET is large compared to Rq. To 
take into account the frequency dependence of the transition rates, includ- 
ing the energy exchange with the measured system, is the main objective 
of the quantum theory presented in section 3.2. 



2.1. TRANSITION RATES 

The dynamics of the SET consists of stochastic transitions between the 
charge states, i.e. by electrons randomly tunneling on and off the island. In 
orthodox SET theory the rates for the different transitions are given by the 
Golden Rule rates for tunneling through the left and right tunnel junctions 

T L / R (E if ) = ^y^^/^, (2) 

where Eif is the energy difference between the initial and final state and 
the Bose function appears from the convolution of the two Fermi functions 
for the filled initial states and the empty final states in the leads. E^ has 
two terms; one is the change in charging energy, and the other is the work 
done by the voltage bias. 

To be specific we now consider the electrons to gain energy from the 

J / Ft 

voltage source by tunneling from left to right. Then the rates T n ± for 
transitions from charge state n to n ± 1 across the left /right junction are 

It + = T L (E n{n+1) + eV L ), T R + = T R (E n{n+1) - eV R ), 

r£_ = r L ( J E ra(ri _ 1) - eV L ), r R _ = T R (E n{n _ 1} + e v R ), { ) 

where E mn = E^m) — E^n) and V L / R = a L / R V is the voltage drop over 
the L/R junction. We denote the sum of rates taking the SET island from 
ntonil with S n± = T^ ± + T R ± . 



2.2. MASTER EQUATION AND STEADY-STATE 

During a specific measurement, i.e. a realisation of the stochastic process, 
the number of (extra) electrons on the SET island as a function of time 
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Figure 2. In orthodox theory the transitions between charge states in the SET are 
considered to be instantaneous. For a specific realisation (measurement) this results in: 

a) The currents through the left/right junctions Il/Ir consist of a train of delta-peaks. 

b) The charge on the island N(t) is a telegraph signal. In order to calculate the transition 
rates one uses the Golden Rule approximation, which includes integration over a timescale 
of h/Eif, where Eif is the energy gain in the transition. This is indicated as grey areas 
in figures a) and b). c) In the real-time Keldysh formalism this timescale of tunneling 
reappears as a timescale on which the density matrix is off-diagonal. 



is piecewise constant, with steplike changes due to instantaneous tunnel 
events (see Fig. 2a-b). Mathematically one may describe the dynamics of 
this process averaged over a large number of measurements. Introducing 
the probability P n (t) of finding the SET in charge state n at time t we may 
write down a master equation 

d t P(t) = ± • P(t), (4) 

where P = [. . .P_i(i) Po(t) Pi(i) . . .] T is a column vector and X is the 
tridiagonal transition rate matrix with off-diagonal elements (S)( n±1 ) n = 
E„±, and to ensure probability conservation (S)„„ = — (S n + + £„_). The 
solution to Eq. (4) may be written P(t) = II(t)-P(0), introducing the time- 
evolution operator (matrix) IL(t) = e s *. To compute the matrix exponent 
one needs to diagonalize S. The steady-state corresponds to the eigenvector 
of S with eigenvalue zero, and the other eigenvalues determine rates for 
exponential decay of the corresponding eigenvector. 



2.3. CHARGE, VOLTAGE AND CURRENT 

In order to calculate physical quantities we introduce operators for the 
number of (excess) electrons on the island (N) and the tunnel currents 
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Figure 3. When an electron tunnels from/to the leads to/from the island the charges 
on the capacitances will redistribute. The charge transport across the boundary of a 
region containing the SET island plus its nearby capacitances is shown for a) an electron 
tunneling onto the island from the left lead, and b) an electron tunneling from the island 
into the right lead. ea L / R represents a displacement charge. 



across the left (II) and right (Ir) junctions, 

(N) nn = n, (I L ) {n±1)n = ± e r£ ± , (I R ) {n±1)n = =Fer^±, (5) 

noting only the non-zero elements, i.e. N is diagonal, and Il/r are both 
tridiagonal with zeros on the main diagonal. The steady-state expectation 
value O st of an operator O is given by 

O st = (6) st = [1...1] -6-P st , P st = ft(oo) • P(0) => (6) 
N st = ^nP n st ,7£* = e^P n ^ + -rt),^ = e^P n s *(r«_-r« + ), 

n n n 

which defines the meaning of the brackets ( ) s t- One gets the voltage of 
the SET island by multiplying the number of extra electrons N by e/Cs, 
and adding the gate bias dependent constant according to Eq. (1). One 
should also note that each tunnel event in the SET is followed by a fast 
redistribution of charge on the nearby capacitances, see Fig. 3. Thus a 
tunnel event in one junction will create a displacement current also in 
the opposite lead. The operator I for the externally measurable current 
is given by[20, 21, 22, 23, 24] I = a L I L + a R I R , (see Eq. (1)). For the 
steady-state current we have I st = Iff = I R , using the detailed balance [25] 
P**S n+ = Pn+i^ln+i)--) but for finite frequency properties the difference 
between tunnel currents and the externally measurable current becomes 
important. 



3. Noise - Fluctuations and Correlations 

The SET dynamics is noisy due to its stochastic nature. The fluctuations in 
the current through the SET determine the measurement time (t m ) needed 
to separate the dc currents corresponding to different charges at the input 
of the SET. The fluctuations of the charge on the SET island induce a 
fluctuating voltage on the capacitance coupling to the measured system 
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(see Fig. 1), which may be disturbed. This effect is called the back-action 
of the measurement, i.e. the meter acting back on the measured system. 

The fluctuations in the SET can be thought of in terms of two contri- 
butions, one from the shot noise of the sequential tunneling described by 
orthodox theory, and one from the quantum fluctuations. At low frequency, 
the shot noise will dominate, as described by orthodox SET theory in sec- 
tion 3.1 below. At high frequencies the SET bias may be neglected and 
the quantum fluctuations, i.e. Nyquist noise, will dominate. The Nyquist 
voltage noise is given by the impedance of the SET island to ground [26], 
i.e. from the two junctions in parallell (Z p (uj)), while the current noise is 
given by the impedance through the SET, i.e. from the two junctions in 
series (Z s (u>)): 



S vv (lo) 



2huRe{Z p (uj)} 
Re{Z s (u;)} 



e 2 R K 
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R 



+ 0(1), 



(7) 



given to leading order in . 

To obtain the spectrum of fluctuations for intermediate frequencies it 
is necessary to solve the full quantum problem, which was done in [27] for 
the voltage noise. The result of this calculation and the comparison to the 
shot noise and quantum fluctuation results in the two limits, are shown in 
Fig. 4 for sample # 1 described in Table 1. As can be seen, the result of the 
full calculation coincides with the shot noise and the quantum noise in the 
low and high frequency limits respectively. 

In addition to this spectrum, the RF excitation gives a component at 
/rf, and the non-linearity of the IV-characteristics gives an additional com- 
ponent at 3/rf- However these frequencies are much lower than the relevant 
mixing frequency AE/h for the SCB-qubit described in section 5. 



3.1. NOISE IN THE SET: ORTHODOX THEORY 

Since orthodox SET theory gives the correct low-frequency limits for the 
noise, and is also the natural reference for discussing the quantum theory, 
we will describe this in some detail. Together with section 2 the following 
discussion is sufficient for calculating the orthodox theory noise, including 
an arbitrary number of charge states. In section 4 we will go through the 
analytically solvable case of two charge states. 

We are interested in the fluctuation of charge and current around their 
average values. Therefore we define fluctuation operators, 50 = O — O st l, 
where 1 is the unit matrix with the same dimension as O and O st is the 
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Figure 4- Comparison of the calculations of the symmetrized spectral density of voltage 
fluctuations of the SET island using the full quantum mechanical calculation (full line), 
the classical shot noise (dashed line), and the quantum fluctuations assuming a linear 
SET impedance (dotted line). Parameters for sample #1 described in Table 1, are used 
in all calculations 



steady-state expectation value of the operator. We may then write the 
following expressions for the fluctuation correlation function (r > 0) 

K AB (T) = (SA-U(r)-SB) st , K ab (-t)=K ba (t). (8) 

The master equation (4) is dissipative and therefore irreversible in time, 
which explains the need for the special negative r definition. Special care 
has to be taken for the autocorrelation of a current pulse, i.e. the shot- 
noise. Since the first current operator instantaneously changes the state of 
the system, the second current operator does not operate on the same state, 
and thus the shot-noise is not included. Using some arbitrary representation 
for the S- function current pulse on the tunneling time-scale Eif/h one gets 
/ 5(t — t )5(t)dt = 5(t ) which gives 

K IrIr {r) - K IrIr (r) + e<5(r) £ P^(T r n+ + F r n _), r G {L, R}. (9) 

n 

We will also need the spectral densities of the fluctuations, i.e. the Fourier 



SET BACKACTION 



9 



transforms of the correlation functions 

/oo rco 
e-^ T K AB {r)dT = / e-^ r K AB {r) + e lu > T K BA (r)dT. 
-oo JO 

(10) 

By e.g. diagonalizing £ one finds that the needed Laplace transform of II(t) 
is given by Eq. (13), just replacing S(cj) — ► S. 

Here we use the unsymmetrized definitions of the correlation functions 
in order to separate noise where the tunnel processes absorb the energy 
\huj\ (positive frequencies) from noise where the tunnel processes emit the 
energy \%uj\ (negative frequencies), when we go to the quantum expres- 
sions in the next section. When this separation is not needed one may use 
the symmetrized definition 5^ l (o;) = Jf^ e~ luJT (Ka B (t) + Kba(t))(It = 
Sab(u) + Sab(— oj). In the classical calculations Sh(uj) = Su(-uj); there- 
fore the symmetric definition only gives a factor of 2 needed to recover 
the usual expression for the shot-noise across a single junction Sfj m (uj) = 
2e/ s *[28]. 

Apart from using unsymmetrised correlation functions, the formalism 
for the orthodox theory presented in this section is equivalent to, and in- 
spired by, the work of Korotkov[22] . 

3.2. NOISE IN THE SET: QUANTUM THEORY RESULTS 

The master equation (4) contains only dissipative transitions between SET 
charge states. The fact that the SET might be in a coherent state like 
\N = 0) + \N = 1) was only taken into account in deriving the tunnel- 
ing rates. This Golden Rule derivation includes an integration over the 
timescale toR = h/Eif, and the main approximation behind Eq. (4) is that 
this timescale is short compared to the dynamics you describe. For steady- 
state or low-frequency properties this is fulfilled since the tunneling rates 
include the small tunnel conductance, i.e. 1/E ~ tcaRx / Rt *C ten- The 
shot-noise term had to be handled separately since it corresponds to corre- 
lations on the time-scale tcR- For high frequency dynamics, i.e. Tiu ~ Eif, 
the approximation breaks down and one has to consider the effects of quan- 
tum coherence. 

One way to include quantum coherence is to use the real-time dia- 
grammatic Keldysh approach described by Schoeller and Sch6n[29]. In the 
sequential tunneling approximation the low-frequency results coincide with 
orthodox theory [22], and the diagrams used may also be compared with 
the time evolution of the state of the SET island, see Fig. 2c. We will not 
further describe the method here, only state the results in such a way that 
the finite frequency noise, including arbitarily many charge states, may be 
calculated. 
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3.2.1. Inelastic transition rates 

When we take into account that the SET can absorb or emit energy the 
transition rates are modified. First of all the tunnel rates become frequency 
dependent, since the noise energy \hu\ should be created or absorbed in a 
tunnel event. The rates 

r£±M = \v L {E n(n±l) ±eV L + nuj), r* ± H = ^r R (E n{n±1)T eV R + hco), 

(ii) 

are here defined so that the transition is facilitated by positive hco, i.e. the 
tunnel event absorbs energy for positive fko. The definition of the S-matrix 
below Eq. (4) is still valid replacing the rates with the sum of positive and 
negative frequency rates 

s^M=r^M + r^(- w ), s n± (u,) = e^M + (12) 

Note that in the zero-frequency limit the orthodox transition rates are 
recovered, i.e. £ n ±(0) = S n ±. The time-evolution is evaluated by Laplace- 
transformation, and the Laplace transform of the time-evolution operator 
obeys the following Dyson type of equation [31] 



1 



-£(„) 



(13) 



where 1 is the unit matrix with the same dimension as the S-matrix, de- 
termined by the number of relevant charge states. 

3.2.2. Finite Frequency Noise Expressions 

We now present the noise expressions valid at finite frequency. The charge 
noise is given by 



S nn (uj) = 2Re [1 . . . 1] ■ N ■ n(w) • iViH • P 



St 



(14) 



where the charge operator now is tridiagonal and frequency dependent when 
it stands in the first position (Ni), 



(Aq) (n±1)n = TtT n ±(u))/u, (Ni) 
The current noise is given by 



-[(^l)(n-l)n + (^l)(n+l)n]- (15) 



5j/H = 5|J ot H + [i...i]-/ 2 H- fl(w) + n(-w) -/iH-p 5 *, (16) 



where the shot-noise contribution is 

sfr t (co) = e 2 £ 2P„VK + M + r;„H], 

n,r£{L,R} 



(17) 
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Figure 5. Calculated Svv{w) (bold solid line) and Sh(lu) (bold dashed line) for a 
symmetric SET at zero temperature with = R§ = Rk and Ec = 3fcs. The thin 
lines are the results from orthodox theory. The SET is biased at n x — 0.25, and slightly 
above the Coulomb threshold eV = 1.1E C , with I st w 0.5nA. Both S V v(u) and Sh(lo) 
go to zero at Tito = 1.05Ec which correspond to the maximal extractable energy. 



and the current operators are both tridiagonal matrices with non-zero en- 
tries 



(lLl)(n±l)n 
{Ihljnn 
(lL2)(n±l)n 



±e 



r£ ± (o) + r£ ± M 



e [r^ + (o)-r^ + H] 
±e^ ± {uj) , (i L2 ) n . 



= 0. 



(18) 



The operators for the right junction are constructed by changing L — > R and 
multiplying by —1. The total current operator is still I n = ail in + (XrIru- 
We see that also the current and charge operators aquire a frequency 
dependence and that the first and second operator have different form, due 
to correct ordering along the Keldysh contour. 



3.3. DIFFERENT LIMITS OF THE QUANTUM EXPRESSIONS 

An example of Svv(v) an d Sji(u) is shown in Fig. 5, with the orthodox 
results as comparison. In the limit to — > the orthodox results are recov- 
ered. The quantum corrections in the low-frequency regime are discussed 
in detail for the two-level system in section 4.2. The noise at large negative 
frequencies, which correspond to noise where the SET emits large ener- 
gies, tend to zero. In our approximation the voltage noise is analytically 
zero when P^Y n ±{u) = for all n, i.e. when no inelastic tunneling events 
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are allowed from the steady-state. The shot-noise part of the current noise 
behaves similarly, and the correction has the same order of magnitude as 
the cotunneling shot-noise, which we do not take into account. Therefore, 
to the accuracy of our approximation, the current noise vanishes when the 
shot-noise vanishes, i.e. when P^T n ±(uo) = for all n. 

In the high positive frequency limit the spectral noise density of the 
SET should be independent of the bias and be dominated by the Nyquist 
noise in Eq. (7). In this limit all rates entering the quantum noise formulas 
simplify to 

r&*M - rv* s J* £ + oil) = + o(i), (19) 

where 0(1) is a bias-dependent constant. Using Syy{u) = e 2 /C|; Snn{uj) 
we thus arrive at the following high-frequency quantum noise expressions 

SwH « 4^ 2 ( rL + rR )' « « 4e 2 (alr L + a 2 R T R ), 

(20) 

where we used that J2 n = 1- One may easily check that this agrees with 
Eq. (7). 



4. Two level approximation for the SET 

As an explicit example we now show the analytically solvable case with low 
bias and temperature, such that only the two lowest energy charge states, 
say N = and N = 1, are occupied, i.e. the only non-zero rates are £o+ 
and Si_. 



4.1. ORTHODOX THEORY 

The master equation Eq. (4) then simplifies to 

with the solution given by the time-evolution operator 

/ pst pst \ / pst pst \ 

n(t) = ( P \t pst) + [ _ l P st ) e ~ Et it > 0), (22) 

where we defined the sum rate £ = £o+ + and the steady-state occu- 
pation probabilities P(f = Ei_/E and Pf = S 0+ /S. The first term in the 
time evolution matrix (operator) Tl(t) gives the steady-state solution, and 
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the second term shows simple exponential relaxation with a single rate S 
towards the steady state. The charge and tunnel current operators are 




1 





1 



1 




(23) 



The steady-state properties and correlation functions are 

N st = pst^ pt = jst = pst^ Q+ = p** Sl _ ) 



K nn (t) = PfPfe'^, S NN ((d) 
Kn(r) = 



21 st P S %+ + PfEi- 



(d 



2 + £ 2 



2 + £ 2 



-St 



= eI st [(a R + ai)5{r) + (a L S 0+ - a R £ 1 _)(a i? £ + - a L £i_) — 
S n ((d) = eP\al + a\) + 2 ( " lS ° + = °fl s i-)("fl s o + - a L Ex. 



E 2 + w 2 



s ^o+ + st + ^ 2 K + "l) 

S 2 + w 2 



(24) 



Notice that both charge- and current-noise are proportional to the steady- 
state current. 



4.2. LOW FREQUENCY QUANTUM THEORY 

In the finite but low-frequency regime, i.e. when still only Tq + (lj) and 
((d) are non-zero, the quantum expressions simplify to: 



~>NN 



((d) 



P<f2T£ + H + Pf*2rf_( 



id 



w 2 + (S 0+ H + S 1 _(o;)) 2 



s n ( 



2e 2 



M) "L 1 0+ + "M a R L 1 



+ 



u; 2 + (£ + + £i-) 2 

e/ st (S 0+ - Ei_) + 2e 2 (S + + Ei_)(PfWf_ - P<fa L r, L 



(25) 



where the frequency dependence of Y, n ±((d) and T^^ R (o;) has been sup- 
pressed for brevity. At zero temperature, with symmetric junctions Cl = 
Cr and Rj> = i?^ = Rt, the expressions simplify further: 



Snn(u) 



Sii(u) = - 



I st + 



I st + 



Hid 
Hid 



1 



2eP T 



(So+W + S^O^+u, 2 ' 

(£ 0+ (0)-£i-(0)) 2 



1 + 



(E 0+ (0) + S 1 _(0)) 2 + ^ 2 



• (26) 
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We find that the classically derived expressions (Eq. 24), symmetric in ±u, 
get an asymmetric quantum correction linear in uj. The correction is added 
to the steady-state current and is proportional to hco/eRx, which would be 
the current through one junction with voltage bias V = hu/e. 

Thus the important difference between the orthodox result and the full 
quantum result is that we get an asymmetry in the noise spectrum be- 
tween the positive and negative frequencies in the full quantum result. 
This becomes important for example in a qubit measurement since it will 
drastically change the occupation of the two levels in the qubit. In the next 
sections we will discuss such a qubit, and how the noise from the SET will 
affect the measurement of the qubit. 



5. The Qubit 

The qubit we consider is made up of the two lowest lying energy levels in 
a single Cooper-pair box (SCB) [17]. An SCB is a small superconducting 
island, with charging energy Eq = e 2 /2(C q b + CJ* + C c ), coupled to a 
superconducting reservoir via a Josephson junction with Josephson energy 
Ej. In order to have a good qubit the following inequalities have to be 
fulfilled: A s > Eq » Ej S> fc^T, where A s is the superconducting energy 
gap and T is the temperature. The low temperature is required to prevent 
thermal excitations and the high superconducting gap is needed to suppress 
quasiparticle tunneling. For suitable values of the gate voltage (close to 
n g = 1/2) the box can be described by the following two-level Hamiltonian 
[17, 15] 

Hf = -^(1 - 2n g )a z - ^-a x (27) 

written in the charge basis (|| = (n = 0| = (1 0), (|| = (n = 1| = (0 1), 
where n is the number of extra Cooper-pairs on the island, a XtZ are the 
Pauli matrices, and n g = C g V g /2e is the number of gate-induced Cooper- 
pairs. By changing the gate voltage the eigenstates of the qubit can be 
tuned from being almost pure charge states to a superposition of charge 
states. The eigenstates of the system written in the charge basis are 

|0> = cosfa/2)|T> + sinfa/2)|i> 

|1> = -sinfa/2)|T> + cosfa/2)|I>, (28) 
where rj = arctan(Ej / 4:E q b(l — 2n g )) is the mixing angle. The energy dif- 
ference between the two states is AE = ^(4E qb ) 2 (l-2n q ) 2 + E 2 j and the 
average charge of the eigenstates is 

Qo = 2e(0| |0) = 2esin 2 (r ? /2), Q x = 2e(l| [){[ |1> = 2e cos 2 (r ? /2). 

(29) 
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6. SET Charge Measurement of Qubit 

We assume that the qubit is in the state co(0)|0)+ci(0)|l) before a measure- 
ment. A perfect charge measurment, i.e. qubit read-out, will now give the 
charge Qq or Qi with probability | Co | 2 and |ci| 2 respectively. The two qubit 
states correspond to slightly different SET gate voltages, and therefore to 
two slightly different sets of transition rates, and two different steady-state 
currents iff and if. The current fluctuates so there is a finite measurement 
time t m needed to separate Iff from If* [15] 

.theory _ i 1 ! ~ ^ f _ e 2 R T Cy, , , 

m ~ 85//(0) ~ AQ 2 ' (6U) 

where AQ = (Qi—Qo)-C c /C q b is the charge difference seen by the SET, and 
where for the last estimate have we used a symmetric SET biased slightly 
above the Coulomb threshold at n x = 0.25. The experimentally determined 
charge sensitivity 5Q gives the measurement time according to 

and by comparing Eq. (30) and Eq. (31) we may deduce a rough theoretical 
estimate for the charge sensitivity 5Q theory . For the parameters of sample 
# 1 and # 2 listed in Table 1 we find 5Q theor y = e^R T C^/2 ~ 1.5 • 
10~ e e/y/ Hz, which further substantiates that the measurements are indeed 
almost shot noise limited. 



7. SET Back Action on the Qubit 

When performing a measurement on a qubit, the measurement necessarily 
dephases the qubit. However, there can also be transitions between the two 
qubit states which destroys the information that the read-out system tries 
to measure. This mixing occurs due to charge fluctuations of the SET island 
Snn(oj), which creates a fluctuating charge on the coupling capacitance C c , 
which is equivalent to a fluctuating qubit gate charge n g — ► n g + Sn g (t). 
The characteristic time for this mixing process is the time t m i x . In the weak 
coupling limit C c <C {Cs, C q b}, which is relevant for the qubit measurement 
setup, we may use pertubation theory to evaluate t m i x . 

The fluctuating term in the qubit Hamiltonian, written in the charge 
basis is 

ATP Jp^C 

5Hf(t) = -^±25n 9 {t)a z = Ej5n(t)a z , Ej = (32) 

where e Sn(t) represents the fluctuations of charge on the SET-island, and 
where we have omitted the term quadratic in 8n g (t). We also defined an 
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electrostatic SET-qubit interaction energy Ej. In the eigenbasis, the full 
qubit Hamiltonian now reads 

AE 

5H e {t) = —a z + E I Sn(t)(cos(rj)a z + sm(r])a x ). (33) 

The cos (77) term leads to dephasing, and the sin(?7) term to interlevel tran- 
sitions (level mixing and relaxation). 

Weak coupling means E T < {E c ,Eg} and we may use standard sec- 
ond order time-dependent pertubation theory [30] to see the effect of the 
charge fluctuations. Assuming that the qubit is in a pure coherent state 
co(0)|0) + ci (0)| 1) at time t = 0, we may express the effect of the fluctu- 
ating SET charge in terms of the time-evolution of the quantities PQ b (t) = 
\c {t)\ 2 ,Pf(t) = |ci(t)| 2 and = |co(*) c * (*) I- In terms of the q ubit 
density-matrix, Pg^(i) are the diagonal elements, determining the occu- 
pation of respective state, and £(t) is the magnitude of the off-diagonal 
elements, describing the quantum coherence between the states. We find 
that the qubit occupation probabilities obey a similar master equation as 
the two-level SET (Eq. (21)), where now the transisition rates are deter- 
mined by the spectral density of the charge fluctuations at the frequency 
corresponding to the qubit level splitting: 

Eft. = ^f^S NN (AE/h), Sf_ = ^h^ls NN (-AE/h), (34) 

where Snn(lu) is the asymmetric charge (number) noise spectral density de- 
fined in Eq. (10). The qubit therefore relaxes exponentially on the timescale 
t mix , where 

tmix = S g6 = E 0+ + (35) 

towards the steady-state P{j b ' st = E^/£« 6 , Pf M = sft/E" 6 . One may 
note that the quantum corrections in Eq. (26) cancel in the expression 
Snn(w) + Snn(-u) determining the mixing time. 

Due to the charge fluctuations the quantum coherence decays exponen- 
tially 

t{t)=me-* /T ", T- 1 = 2^±j^S NN (0) + ^, (36) 
where t v is the timescale for dephasing. 

7.1. SIGNAL TO NOISE RATIO IN QUBIT READ OUT 

Now we can use the measured data for samples #1 and #2 to calculate both 
t m and t m i x , and thus we can also get the expected signal-to- noise ratio for 
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Figure 6. Calculated Sy™(f) for samples #l(full line) and #2(dashed line). The arrows 
show the energy separations of the two qubit states for an aluminum and a niobium qubit, 
respectively. 



a single-shot measurement, which is simply given by SNRss = ytmkjtm- 
The spectral densities Svv(f) for the two samples at the optimum charge 
sensitivity, at a current of 6.7 and 8.0 nA for samples #1 and #2, re- 
spectively, are displayed in Fig. 6. As can be seen both through Eq. (34) 
and in Fig. 6, the mixing time increases strongly with increasing AE. In 
our case AE is limited by the superconducting energy gap A, which for 
aluminum films corresponds to about 2.5 K. Using niobium as the qubit 
material would substantially increase the mixing time. If we assume that 
Eq and thus also AE can be scaled with A and that the coupling C c /C q b 
is kept constant, mixing times of several ms can be reached. In that case, 
other sources than the SET noise would most probably dominate the mix- 
ing. The results for the samples #1 and #2 are summarized in Table 1, 
where a coupling C c /C q b = 0.01 is assumed. 

7.2. COULOMB STAIRCASE 

One may use the SET to measure the so-called Coulomb staircase, i.e. the 
average charge of the qubit as a function of gate voltage. In an ideal situa- 
tion with no energy availible from an external source, at zero temperature, 
the qubit would follow the ground state adiabatically and the charge would 
increment in steps of 2e at n g = n + 0.5, n integer. These steps are not per- 
fectly sharp because of the Josephson energy mixing the charge states[17]. 
Now assuming that the qubit equilibrium is determined by the SET back 
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Figure 7. The calculated Coulomb staircase of an SCB as measured by the SET, as- 
suming that the SET backaction determine the steady-state of the qubit. 



action the charge measured should instead be Pq ' s Qo+Pf ' s Q\. The steps 
are now rounded further [32] due to the finite probability for the qubit to be 
excited, see Fig. 7. The quantum asymmetry of the noise-spectrum, indicat- 
ing the difference of qubit excitation and relaxation, is vital to recover the 
correct Coulomb staircase. Straighforward use of orthodox theory would 
predict an equal population of the qubit states. 
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TABLE 1. 



SET 


Qubit 


SQ 


AE 




Cm 


tmix 


SNR 


Sample 


material 


e/v^S 


[K] 


[nV 2 /Hz] 


M 


M 




1 


Al 


6.3 


2.4 


0.29 


0.40 


8.6 


4.6 


1 


Nb 


6.3 


15.5 


0.056 


0.40 


1860 


68 


2 


Al 


3.2 


2.4 


0.39 


0.10 


6.4 


8.0 


2 


Nb 


3.2 


15.5 


0.080 


0.10 


1300 


114 
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